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g^: ABSTRACT 

We present results of a simulation of two flavor QCD on a 16^ x 4 lattice using p4-improved 
^ ' staggered fermions with bare quark mass m/T = 0.4. Derivatives of the thermodynamic 
grand canonical partition function Z(V, T, /i^j, /i^) with respect to chemical potentials iJiu,d 
for different quark flavors are calculated up to sixth order, enabling estimates of the 
pressure and the quark number density as well as the chiral condensate and various sus- 
ceptibilities as functions of yUg = (/iu + /irf)/2 via Taylor series expansion. Furthermore, we 
analyze baryon as well as isospin fluctuations and discuss the relation between the radius 
of convergence of the Taylor series and the chiral critical point in the QCD phase dia- 
gram. We argue that bulk thermodynamic observables do not, at present, provide direct 
evidence for the existence of a chiral critical point in the QCD phase diagram. Results 
are compared to high temperature perturbation theory as well as a hadron resonance gas 
model. 



PACS numbers: 12.38. Gc, 12.38.Mh 



1 Introduction 



The thermodynamics of strongly interacting matter has been studied extensively in lattice 
calculations at vanishing quark chemical potential, /ig (or baryon /ib = '^l^'q) [!]• Current 
lattice calculations strongly suggest that the transition from the hadronic low temperature 
phase to the high temperature phase is a continuous, non-singular but rapid transition 
happening in a narrow temperature interval around the transition temperature Tq ~ 
170 MeV. Recent advances in the development of techniques for lattice calculations at 
non-zero quark chemical potential /ig [2-6] have also enabled the first exploratory studies 
of the QCD phase diagram and bulk thermodynamics in a regime of small /x^, i.e. for 
^lJT < 1 and T > 0.8 %. 

Guided by phenomenological models which suggest that at low temperature and non- 
zero quark chemical potential the low and high density regions will be separated by a first 
order phase transition it has been speculated [7] that a 2'"^ order phase transition point, 
the chiral critical point, exists in the interior of the QCD phase diagram, at which the line 
of first order transitions ends. For smaller values of ^J^q/T the low and high temperature 
regime will then be only separated by a crossover transition. The first exploratory stud- 
ies at non-zero quark chemical potential indeed gave evidence for such a chiral critical 
point [2], although subsequent investigations made clear that at present any quantitative 
statement about the location [8, 9] and maybe even about the existence of such a chiral 
critical point is premature. The first calculations of the baryonic contribution to the 
pressure in strongly interacting matter [3, 10-14] also suggested that the transition and 
the thermodynamic behaviour in the high-T phase resemble the picture found previously 
at zero net baryon density. In this case (/i^ = 0) the transition occurs over a narrow 
temperature interval. Beyond a region T e [To, I.STq] where there are large deviations 
from ideal gas behaviour, thermodynamic observables rapidly approach the high-T ideal 
gas limit; for instance, the pressure agrees with the Stefan-Boltzmann prediction to within 
~ 20%. In the low temperature hadronic phase it has been found that a hadron reso- 
nance gas model provides an astonishingly good description of basic features of the T and 
//g-dependence of thermodynamic observables [15]. 

Maybe one of the largest changes compared to the thermodynamics at //g = has 
been observed in the temperature dependence of the quark number and isovector sus- 
ceptibilities [11]. For = these observables show a similar temperature dependence. 
They have been found to change rapidly at the transition temperature but continue to 
increase monotonically at larger temperatures [16-20]. For Hq > 0, however, the quark 
number susceptibility develops a pronounced peak at the transition temperature while 
the isovector susceptibility continues to show a temperature dependence similar to that 
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found at /x^ = 0. Such a behavior, indeed, is expected to occur in QCD in the vicinity 
of a 2"*^ order phase transition point [21]. Susceptibihties thus may provide the most 
direct evidence for the existence of a 2"*^ order phase transition in the QCD phase dia- 
gram. Finding the characteristic volume and/or quark mass dependent universal scaling 
behavior*^ of susceptibilities would undoubtedly establish the existence of a chiral critical 
point. 

In this paper we want to extend our previous study of thermodynamics at non-zero 
quark chemical potential [3,11], which is based on a Taylor series expansion around fig = 0, 
to the 6*^ order^. Going to higher orders in the expansion is of particular importance 
for the analysis of higher derivatives of the density of the grand potential expressed in 
units of the temperature^, Q(y,T, /x) = {VT^)~^lnZ, i.e. for an analysis of generalized 
susceptibilities. Accordingly our main emphasis will be to further analyze the properties 
of quark number and isovector susceptibilities, relate them to diagonal and non-diagonal 
flavor susceptibilities, calculate the chiral susceptibility and discuss to what extent the 
pronounced peaks found in some of these susceptibilities give evidence for the existence 
of a chiral critical point in the QCD phase diagram. 

This paper is organized as follows. We start in the next section by summarizing 
basic results on QCD thermodynamics at non-zero quark chemical potential obtained in 
high temperature perturbation theory [22, 23] and expectations based on properties of 
the hadron resonance gas model at low temperature [24,25]. In section 3 we present 
results on the calculation of various thermodynamic observables obtained from a Taylor 
expansion up to G*'^ order in jJLq/T. In section 4 we analyze the convergence properties 
of the Taylor series. Section 5 is devoted to a discussion of the reweighting approach to 
QCD thermodynamics at /Xg 7^ and its comparison to the Taylor expansion approach. In 
section 6 we give our conclusions. An Appendix contains details on the Taylor expansion 
of various observables studied here. 



2 Thermodynamics at low and high temperature 

In this section we want to briefly discuss the density or chemical potential dependence of 
thermodynamic observables in the asymptotic high temperature regime as well as at low 

temperatures in the hadronic phase of QCD. On the former we have information from 
high temperature perturbation theory which recently has been extended to 0{g^ In 5') [23] 
also for non-vanishing quark chemical potentials /xj, f = u, d, .., and as for Hf = [26] 

''The 2"'' order critical point in the QCD phase diagram is expected to belong to the universality class 
of the 3-dimensional Ising model. 

''Some results on the radius of convergence have been reported recently in an 8*'' order Taylor expansion 
for two flavor QCD [20]. 

■^We will call this in the following the grand potential although conventionally the extensive quantity 
VT'^Cl is called the grand potential. 
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is thus now known to all perturbatively calculable orders. Although a comparison of 
the perturbative expansion with lattice calculations at yUj = suggests that quantitative 
agreement cannot be expected at temperatures close to the transition temperature, Tq, 
we can gain useful insight into the structure of the Taylor expansion used in lattice 
calculations to study thermodynamics at /i/ 7^ 0. 

In the low temperature hadronic phase a systematic QCD based analysis is difficult 
and one generally has to rely on calculations within the framework of effective low energy 
theories or to phenomenological approaches. Here we will focus on a discussion of the 
properties of a hadron resonance gas model, which recently has been compared to lattice 
results at non-zero quark chemical potential quite successfully [15] and which also is known 
to describe experimental results on the chemical freeze-out of particle ratios observed in 
heavy ion collisions rather well [25]. 



2.1 High temperature perturbation theory 

In the infinite temperature limit the grand potential of QCD*^, Q{T, /i) = In Z{V, T, fi) /VT^, 
which is equivalent to the pressure in units of T^, approaches that of a free quark-gluon 
gas (Stefan-Boltzmann (SB) gas). 



PSB 
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(2.1) 



where the first term gives the contribution of the gluon sector and the sum over the 
fermion sector extends over Uf different flavors. In Eq. (|2.H) we only gave the result 
for massless quarks and gluons. Also in the following we will restrict our discussion of 
perturbative results to the case of QCD with massless quarks. In this section we also 
use 11 to denote the entire set of Uf different chemical potentials (/i„, fid,---)- We also 
introduce the shorthand notation /x^ = X)/ (z^//^)^- 

The additive structure of contributions arising from gluons and the different fermion 
flavor sectors persists at 0{g'^). Only starting at 0{g^) is a coupling among the different 
partonic sectors introduced. At this order it is induced through the non- vanishing electric 
(Debye) mass term niE [23], 



n{T, /i) = w (T, /i) + / ^) ^ ^3 ^(3) ^rp^ ^) ^ (^(^4) ^ ^2.2) 

'^We suppress here the volume dependence of Vl. Perturbative calculations are performed in the 
thermodynamic limit. The volume dependence, however, has to be analyzed more carefully in lattice 
calculations. 
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with 
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f=u,d,. 

The electric mass term introduces a dependence of a given quark flavor sector on changes 
in another sector, i. e. the quark number density in a flavor sector i, 

TH ^ dn{T, fi) 

T3 dfie/T ' ^ ' ' 

depends on the other quark chemical potentials only at 0{g^). This is also reflected in 
the structure of diagonal and non-diagonal susceptibilities [4, 16-20]. 

d{i,f/Tr ' d{i,f/T)dMT) ■ ^ ■ ' 

The diagonal susceptibilities are non-zero in the ideal gas limit and, moreover, the leading 
order perturbative term stays non-zero also in the limit of vanishing quark chemical 
potential, 

MSM . l + . (2.6) 

The non-diagonal susceptibilities, however, receive non-zero contributions only at 0{g^). 
The leading perturbative contribution is positive and inversely proportional to the electric 
screening mass. However, it vanishes in the limit of vanishing chemical potentials. In this 
case the flrst non-zero contribution arises at 0{g^lnl/g) [22], 

- . (2.7) 



r2 27r5 V 6 27r2'^ / T T 

. (2.8) 



Xfk{T,Q) _ 5 ^6 



^2 144^^6 



As we are going to discuss lattice calculations at non-zero chemical potential which are 
based on a Taylor expansion of the grand potential in terms of /x//T it also is instructive 
to consider a Taylor expansion of the perturbative series for Q. While perturbative terms 
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up to 0{g'^) only contribute to the series up to higher order terms in the expansion 

start receiving non- vanishing contributions at 0{g^). These again arise from the electric 
mass term, i.e. from an expansion of niE/gT in powers of fi^; 

1 + ^V'^' (/x^f + . (2.9) 



7687r7 V 6 

Note that the expansion coefficients up to and including are positive. Starting 

at 0{fi^) they alternate in sign. From the sign of the /x^-term it follows via Eq. ()2.5j] 
that the leading perturbative contribution to the expansion coefficient of diagonal as well 
as non-diagonal susceptibihties at is negative, with the latter being an order of 

magnitude smaller*^. 

Rather than discussing the thermodynamics of QCD at non vanishing quark chemical 
potential in terms of chemical potentials related to quark numbers in different flavor 
channels, it is convenient to introduce chemical potentials which are related to conserved 
quantum numbers considered at low energy, i.e. quark or baryon number and isospin. 
In the case of two flavor QCD, which we are going to analyze in our lattice calculations, 
we thus introduce also the quark chemical potential, fig = + fid)/'^ and the isovector 
chemical potential^ /i/ = {fiu — fJ'd)/'^- In analogy to Eq. ()2.5|) we then can introduce quark 
number and isovector susceptibilities, 

Xq _ _ o . , ^ Xl _ _ o . _ ^ . ^^ 

~ d{n /T)2 ~ v^"" ' rp2 ~ d{ixi/TY ~ ^'^'^^ ' \'^-^^) 

where in the second equality we have assumed degenerate (m, (i)-quark masses. 



2.2 Low temperature hadron resonance gas 

The success in describing particle abundance ratios observed in heavy ion experiments 
at varying beam energies in terms of equilibrium properties of a hadron resonance gas 
model [25] begs a comparison of this model for the low temperature hadronic phase with 
lattice QCD calculations. Indeed this led to astonishingly good agreement [15]. 

In the hadron resonance gas (HRG) model it is assumed that for T < Tg the QCD 
partition function can be approximated by that of a non-interacting gas of hadron reso- 
nances, either bosonic mesons or fermionic baryons. This, however, does not mean that 

''It has also recently been pointed out in the context of a large- TVc expansion that flavor non-diagonal 
contributions to the free energy are suppressed by 0{1/N^) [27]. 

^This definition of the isovector chemical potential differs from that used in [11] by a factor of 2. 
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interactions in dense hadronic matter have been ignored; in the spirit of Hagedorn's boot- 
strap model [24] the inclusion of heavy resonances as stable particles also takes care of 
the interaction among the hadrons in the dense gas at low temperature. 

The partition function of the hadron resonance gas may be split into mesonic and 
baryonic contributions, 

In Zhrg{T,V, ^iq, fii) = lnZ^^(T, + In V, /Xg, /zj) , (2.11) 

iG mesons baryons 

where 

T/ roo 

lnZ*y^(T,l^/i„/i,) = T— dkein{lTz.e-'^^^) , (2.12) 

with energies ef = k"^ + mf and fugacities 

Zi = exp ((35,/i, + 2hifii)/T) . (2.13) 

Here Bi is the baryon number and I^i denotes the third component of the isospin of the 
species in question. The upper sign in Eq. fl2.12j) refers to bosons, and the lower sign 
to fermions. Note that with this convention anti-particles must be counted separately 
in Eq. ()2.11|) with fugacity z~^, and self-conjugate species such as n^,!]' have z = 1. If 
the logarithms are expanded in powers of fugacity, the integral over momenta, k, can be 
performed. This yields 

where the upper factor in braces applies to bosons and the lower to fermions, and K2 
is a modified Bessel function. For large argument, i.e. for mi ^ T the Bessel function 
can be approximated by K2{x) ~ Jti /2x e~^(l -|- Ih/^x -\- 0{x~'^)). Terms with £ > 2 in 



the series given in Eq. ()2.14|) thus are exponentially suppressed. For temperatures and 
quark chemical potentials less than a typical scale of about 200 MeV it generally suffices 
to keep the first term in the sum appearing in Eq. (|2.14j) . The only species for which this 
step would need further justification is the pion; clearly for realistic pion masses more 
care must be taken when evaluating the sum over £. This, however, does not infiuence 
the /iq-dependence of the hadron resonance gas; i?j = for mesons and their contribution 
thus is independent of Hq. At /ij = one readily derives. 



p{T^_ 1 



G(T) + F(T) cosh ( ^ ) , (2.15) 
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with G{T) = (VT^) ^ J2i£mesons^^ F(T) giveii by the Boltzmann approximation 
to the fermion partition function of baryons, 

nr) ^ 1 i: {^p<.{^) . (2..) 

i£ baryons 

Note that each term in the sum for F now counts both baryon and anti-baryon. As the 
meson sector of the partition functions is independent of fig the mesonic component does 
not contribute to Xg, 

|| = 9F(r)cosh^ . (2.17) 

Similarly, at /x/ = we find for the isovector susceptibility 

|^ = G^(T) + F^(T)cosh^ (2.18) 



with 



t£ mesons £=1 ^ 

F\T) = {2h,f(^)\j^) . (2.19) 

is baryons 

where again only in the baryon sector ^ > 1 terms have been neglected. 

The resonance gas model in the (partial) Boltzmann approximation given by Eqs. ()2.15|) - 
()2.19|) leads to simple predictions for the dependence of thermodynamic observables on 
the quark chemical potential /i^. In particular, it predicts that ratios of the density 
dependent part of thermodynamic observables are insensitive to details of the hadronic 
mass spectrum. Nor do they depend explicitly on temperature, but instead only on the 
ratio fiq/T. For instance, one finds 

p{T,fiq) - piT.Q) cosh(3/ig/T) - 1 riq T , f^l^q 



T2 Q IT\ ' tanh ^ . (2.20) 

XqT^ 9cOSh(3/ig/T) flqXq 3/ig \T J 

Similarly ratios of Taylor expansion coefficients of thermodynamic quantities, X, are 
temperature and spectrum independent. For an observable X, of generic form X = 
+ F-^ cosh(3/ig/T), the expansion in Hq/T is given by 

oo 

X = Y.^niT){i.q/T)\ (2.21) 

71=0 
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with 




+F- 
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F 



(2.22) 



2n+l ~ 



0. Hence, ratios are given by 
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for n > 1 



(2.23) 




(2n + 2)(2n + l) 



These ratios as well as ratios of physical observables calculated within the resonance gas 
approximation will be compared to corresponding lattice results in the following. 



3 Taylor expansion for two flavor QCD 

The basic concepts of our Taylor expansion approach to QCD thermodynamics at non- 
zero quark chemical potential have been introduced in [11] where observables have been 
analyzed up to 0{fig). Here we extend the analysis up to the sixth order in fJ,q/T with 
significantly improved statistics. 

Our calculations have been performed for two flavor QCD on an 16^ x 4 lattice with 
bare quark mass ma = 0.1 using Symanzik- improved gauge and p4- improved staggered 
fermion actions. These parameters are identical to those used for our analysis of thermo- 
dynamic observables at non-zero chemical potential up to 0{fiq) [3,11]. The simulation 
uses the hybrid R molecular dynamics algorithm, and measurements were performed 
on equilibrated conflgurations separated by 5 units of molecular dynamics time r. The 
gauge couplings, /3 = used for our calculations cover the interval [3.52,4.0], which 

corresponds to a temperature range T/Tq G [0.76,1.98], where Tq is the pseudo critical 
temperature at /i = for which we use^ Pc = 3.65. This value is used to deflne the 
temperature scale T /Tq. The number of conflgurations generated at each /?- value is given 
in Table ITTl The third and seventh columns give the sample sizes used in [11], where 
coefficients up to n = 4 were calculated. The numbers of additional conflgurations gen- 
erated for the present study of expansion coefficients up to = 6 are listed in the fourth 
and eighth columns. It can be seen that we have increased our statistics in the hadronic 
phase by a factor 4-5 and in the plasma phase by a factor 3. 

For the calculation of various operator traces we use the method of noisy estimators. 
We generally found that expectation values involving odd derivatives of IndetM with 

^In [3] we determined as critical coupling [ic = 3.649(2). Our current analysis favors a slightly larger 
value, Pc = 3.655(5). 
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p 


T/T, 


#(2,4) 


#(2,4,6) 


/3 


T/T, 


#(2,4) 


#(2,4,6) 


3.52 


0.76 


1000 


3500 


3.70 


1.11 


800 


2000 


3.55 


0.81 


1000 


3500 


3.72 


1.16 


500 


2000 


3.58 


0.87 


1000 


3500 


3.75 


1.23 


500 


1000 


3.60 


0.90 


1000 


3800 


3.80 


1.36 


500 


1000 


3.63 


0.96 


1000 


3500 


3.85 


1.50 


500 


1000 


3.65 


1.00 


1000 


4000 


3.90 


1.65 


500 


1000 


3.66 


1.02 


1000 


4000 


3.95 


1.81 


500 


1000 


3.68 


1.07 


800 


3600 


4.00 


1.98 


500 


1000 



Table 3.1: Sample size at each (3 value. 

respect to // arc noisier and require averages over more random vectors than needed 
to estimate expectation values involving only even derivatives. Odd derivatives have 
to appear in even numbers in an expectation value in order for this to be non-zero. 
Such expectation values behave very much like susceptibilities and receive their largest 
contributions in the vicinity of Tq. Still their total contribution to the expansion of e.g. 
the pressure is found to be small up to C(/x^). It only becomes sizeable at 0{iJp). For 
these reasons we used 100 stochastic noise vectors to estimate operator traces on each 
configuration for e [3.60,3.68]. For other /? values 50 noise vectors were found to 
suffice. 



3.1 Pressure, quark number density and susceptibilities 

To start the discussion of lattice results on thermodynamics of two flavor QCD for small 
values of the quark chemical potential we will present results obtained from a Taylor 
expansion of the grand potential Q{T, fi^, fiii) = ^{T,fiq + ^u^q — /i/) and some of its 
derivatives. We will consider expansions in terms of fiq/T at fixed, vanishing iij. The 
pressure is then given by 



;^^Q(r,/.„/.,) = Ecn(T)(^ 
n=0 ^ ^ 



(3.1) 



CP symmetry implies that the series is even in /Xg, so the coefficients c„ are non-zero only 
for n even, and are defined as 



Cn{T) 



n\ d{iiq/Ty 



1 Nl 9"lnZ 



IJ-q 



=0 n\ Nl diiiN^y 



(3.2) 



10 



T/T, 


C2 


C4 X 10 


ce X 10^ 


4 


ci X 10 


4 X 10^ 


0.76 


0.0243(19) 


0.238(61) 


-1.12(121) 


0.0649(6) 


0.098(5) 


0.23(9) 


0.81 


0.0450(20) 


0.377(64) 


1.98(141) 


0.0874(8) 


0.140(6) 


0.44(10) 


0.87 


0.0735(23) 


0.506(68) 


1.69(155) 


0.1206(11) 


0.216(8) 


0.60(13) 


0.90 


0.1015(24) 


0.765(72) 


2.06(159) 


0.1551(14) 


0.302(12) 


0.83(18) 


0.96 


0.2160(31) 


1.491(135) 


4.96(260) 


0.2619(21) 


0.564(23) 


1.47(37) 


1.00 


0.3501(32) 


2.133(121) 


-5.00(359) 


0.3822(26) 


0.839(28) 


0.26(49) 


1.02 


0.4228(33) 


2.258(118) 


-4.49(312) 


0.4501(27) 


0.909(28) 


0.02(44) 


1.07 


0.5824(23) 


1.417(62) 


-5.73(158) 


0.5972(21) 


0.741(17) 


-0.75(26) 


1.11 


0.6581(20) 


0.951(39) 


-1.65(62) 


0.6662(18) 


0.618(11) 


-0.18(10) 


1.16 


0.7091(15) 


0.763(24) 


-0.31(26) 


0.7156(14) 


0.564(6) 


-0.03(4) 


1.23 


0.7517(16) 


0.667(23) 


-0.44(23) 


0.7573(13) 


0.527(5) 


-0.06(3) 


1.36 


0.7880(11) 


0.572(12) 


-0.09(11) 


0.7906(9) 


0.495(3) 


-0.03(1) 


1.50 


0.8059(10) 


0.539(10) 


-0.17(7) 


0.8076(7) 


0.477(2) 


-0.05(1) 


1.65 


0.8157(8) 


0.499(7) 


-0.13(8) 


0.8169(7) 


0.461(2) 


-0.05(1) 


1.81 


0.8203(8) 


0.497(7) 


-0.11(6) 


0.8218(6) 


0.452(1) 


-0.05(1) 


1.98 


0.8230(7) 


0.473(6) 


0.03(4) 


0.8250(6) 


0.441(1) 


-0.03(1) 



Table 3.2: Taylor expansion coefficients c„(T) and c^(T). 



where in the second equality we have explicitly specified that the calculations have been 
performed on a A^^ x A^,- lattice with dimensionless quark chemical potential fi = fiqa. 
Z denotes the lattice regularized partition function for two fiavor QCD. Similarly we 
calculate the quark number density 



as well as the quark and isovector susceptibilities using Eq. ()2.10p . 

= 2C2 + 12C4(^)' + 30C,(^)V... (3,4) 
= 24 + 12cl(^)%30ci(^)* + ... (3.5) 



where 



1 VL{T, fig + fig - Hi) 



(3.6) 
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Figure 3.1: The Taylor expansion coefficients c„ and for n = 2,4 and 6 as functions of 
T/To. 



Explicit expressions for C2, C4, and Cg, C4 and Cg are given in the Appendix. Note that 
the expansion for the quark number susceptibility Xq given in Eq. ()3.4|1 is a derivative 
of the grand potential at yU/ = and thus has the same radius of convergence as that 
of the pressure and quark number density given by Eqs. p.lj) and ()3.3p . The expansion 
coefficient C2 also defines the first term in an expansion of the pressure at non-zero isospin. 
This series may have a different radius of convergence [28]; indeed, since the lightest 
particle carrying non-zero isospin J3 in the hadronic phase is the pion, we might expect 
the expansion to break down in the chiral limit for arbitrarily small /i/. 

The coefficient cq{T) gives the pressure in units of at vanishing baryon density and 
can be calculated using the integral method [29]. It is the only expansion coefficient which 
also requires lattice calculations at zero temperature. Higher order terms can be calculated 
directly from gauge field configurations generated on finite temperature lattices. They, 
however, require additional derivatives of In detM, where M is the quark matrix. They are 
evaluated at fixed temperature, i.e. fixed gauge coupling (3, by calculating combinations 
of traces of products of d"^M/dfi''^ and (see Appendix). 

Results for the Taylor expansion coefficients are listed in Table 13.21 In Fig. ()3.1|) we 
plot Cn and for n = 2, 4 and 6 as functions of T. A comparison with Figs. 3 and 8 
of Ref. [11] reveals the improvement in statistics of the current study. The same features 
are apparent: namely C2 and C2 both rise steeply across To with > C2 as is obvious from 
the explicit expressions given for these coefficients in the appendix; they reach a plateau 
at approximately 80% of the value nf/2 predicted in the Stefan-Boltzmann (SB) limit, 
i. e. for free massless quarks; C4 rises steeply to peak at T ~ Tq before approaching its SB 
limit value nf/47r^ from above, whereas the peak in is much less marked^. 



'^The difference is largely due to the dominance of the disconnected term ((9^ lndetil//9/i^)^) — 
{d^ lndetM/9/i^)^ which contributes to C4 with a coefficient three times that of its contribution to C4. 
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(a) (b) 

Figure 3.2: The /ig dependent contribution to the pressure (left) and the quark number 
density (right) as functions of T/Tq for various values of the quark chemical potential 
calculated from a Taylor series in 6*^ order. Also shown as dashed lines are results from 
a 4*^^ order expansion in /ig/T. 



As can be seen in Table 13.21 in the high temperature phase the 6*^ order expansion 
coefficients generally are an order of magnitude smaller than the 4*^ order coefficients. In 
the low temperature phase they are still a factor 3-5 smaller. As a consequence the Q*^ 
order contributions to the pressure and quark number density are small for /i^/T < 1. 
This is seen in Fig. (13. 2|) which shows Ap/T^ = (p(T, /i^) — p(T, 0))/T^ and riq/T^ in 
the range < ^q/T < 1. Here we also show as dashed lines results obtained from 
a Taylor expansion which includes only terms up to 4*'' order in /Xg/T. This suggests 
that the expansion for the pressure and quark number density is converging rapidly for 
/ig/T < 1. Even for /ig/T = 1 the differences between the 4*'* and 6*'*-order results are 
small and partly influenced by statistics. We also note that in the high temperature 
regime, T>1.5To, our results are compatible with the continuum extrapolated (quenched) 
results obtained with an unimproved staggered fermion action [13]. This supports the 
expectation that deviations from the continuum limit are strongly suppressed with our 
improved action. 

Next we turn to a discussion of quark number and isovector susceptibilities which are 
shown in Fig. (j3.3p . They have been obtained using Eqs. (j3.4|) and (j3.5|) . Again we show 
the corresponding 4*''-order results as dashed lines in these flgures. These lines agree 
with our old results shown as Fig. 9 of [11]. The effect of the new term proportional to 
Cq{T) is to shift the apparent maximum in Xq{T), arising from the sharply peaked /ig- 
contribution proportional to C4(T), to lower temperature. This suggests that the transition 
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Figure 3.3: The quark number susceptibility Xq/T'^ (left) and isovector susceptibility 
Xi/T'^ (right) as functions of T/Tq for various [J^q/T ranging from [iq/T = (lowest curve) 
rising in steps of 0.2 to fiq/T = 1, calculated from a Taylor series in 6*^ order. Also shown 
as dashed lines are results from a 4*'' order expansion in /i^/T. 

temperature at non-zero /ig determined from the peak position of susceptibilities indeed 
moves to temperatures smaller than the transition temperature Tq determined at fiq = 0. 
The figure, however, also shows that at least for T < Tq the G^'^-order contribution can 
be sizeable and still suffers from statistical errors. Better statistics and the contribution 
from higher orders in the Taylor expansion thus will be needed to get good quantitative 
results for susceptibilities in the hadronic phase. 

There is also a pronounced dip in Xq{T) for T/Tq ~ 1.05 which, together with the 
increased error bars makes the presence of a peak in Xq less convincing then it is without 
the inclusion of the /Xg-contribution. However, the error bars also reflect the problem we 
have at present in determining this additional contribution with sufficient accuracy to 
include it in the calculation of higher order derivatives of the partition function. On the 
other side. Fig. ()3.3j) confirms that a significant peak is not present in the isovector channel. 
In fact, if a critical endpoint exists in the (T, /x)-plane of the QCD phase diagram, this 
is expected to belong to the Ising universality class, implying that exactly one 3d scalar 
degree of freedom becomes massless at this point. Since both ijjijj and ipjo'^ are isoscalar 
and Galilean scalars, both are candidates to interpolate this massless field, and hence we 
can expect divergent fiuctuations in both quark number and chiral susceptibilities at this 
point. The latter will be discussed in section 3.3. 

The difference in the temperature dependence of Xq Xi also reflects the strong 
correlation between fluctuations in different flavor components. This will become clear 
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from our discussion of flavor diagonal and non-diagonal susceptibilities in the next section. 

3.2 Flavor diagonal and non-diagonal susceptibilities 

Using the relation between quark number and isovector susceptibilities on the one hand 
and diagonal and non-diagonal susceptibilities on the other hand (Eq. ()2.1()j] ) we also can 
define expansions for the latter, 

2cr+i2cr(^)'+30cr(^)V--- , (3.7) 

2cr + 12crl^^y + 30cfl^^y + --- , (3.8) 

with cr = (c„ + ci)/4 and = (c„ - c^)/4. 

As discussed in the previous section the expansion coefficients c„ and become quite 
similar at high temperature. This was to be expected from the discussion of the structure 
of the high temperature perturbative expansion given in section ITTl as c„ and differ only 
by contributions coming from non-diagonal susceptibilities, which enter with opposite sign 
in these two coefficients. It thus is instructive to analyze directly the expansion coefficients 
of Xuu and Xud- These are listed in Table and plotted in Fig. ()3.4|) . We note that the 
errors on these quantities have been obtained from an independent jackknife analysis and 
thus are not simply obtained by adding errors for c„ and c^. 

Fig. ()3.4|) clearly shows that for T > Tq the various expansion coefficients rapidly 
approach the corresponding ideal gas values, which is zero for all non-diagonal expansion 
coefficients, c^*^. In fact, as discussed in section ITT] the latter receive contributions only 
at Olg^lnl/g"^) for n = 2 and 0{g^) for n > 2. Moreover, it is interesting to note that 
despite the small magnitude of these contributions the leading order perturbative results 
correctly predict the sign of all expansion coefficients for T > To, z.e. cfe < 0, cf > 0, 
and > 0, Cg" < 0. Furthermore, the order of magnitude for 02'^, i.e. {c^'^l — 5 ■ 10~^ at 
T ~ 2Tq, agrees with the perturbative estimate^ [22]. 

A striking feature of the expansion coefficients and c^'^ is that for n > 2 they 
become similar in magnitude close to Tq. In fact, C4" and c^'^ both have pronounced peaks 
at To with (c^'^/c7)peak ^ 0.45 and the expansion coefficients for n = 6 are identical 
within errors. This suggests that any divergent piece in which could occur when 
fiq/T approaches the radius of convergence of the Taylor expansion, will show up with 

'A previous analysis [18] reported values for d^"^ consistent with zero for temperatures T/Tq > 1.25 
within an error of 10^^. This has been found subsequently to be incorrect. The corrected values [20] are 
qualitatively consistent with our findings. (The second panel of Fig. 10 in [20] should be compared to 
twice the value of cf'' shown in our Fig. H3.4|) .l Similar agreement is found with the results of [19]. 
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T/T, 


c^" X 10^ 


cl" X 10^ 


X 10^ 


cf X 10^ 


cf X 10^ 


cf X 10^ 


0.76 


2.23(6) 


0.84(16) 


-0.22(32) 


-1.015(42) 


0.35(14) 


-0.34(29) 


0.81 


3.31(6) 


1.29(17) 


0.60(37) 


-1.060(45) 


0.59(15) 


0.38(33) 


0.87 


4.85(8) 


1.81(18) 


0.57(42) 


-1.177(46) 


0.72(16) 


0.27(36) 


0.90 


6.41(9) 


2.67(20) 


0.72(44) 


-1.339(45) 


1.16(16) 


0.31(36) 


0.96 


11.95(12) 


5.14(39) 


1.61(74) 


-1.148(45) 


2.32(29) 


0.87(57) 


1.00 


18.31(14) 


7.43(37) 


-1.19(101) 


-0.802(32) 


3.23(24) 


-1.31(78) 


1.02 


21.82(15) 


7.92(36) 


-1.12(88) 


-0.681(27) 


3.37(23) 


-1.13(69) 


1.07 


29.49(11) 


5.39(20) 


-1.62(46) 


-0.369(17) 


1.69(11) 


-1.25(33) 


1.11 


33.11(9) 


3.92(12) 


-0.46(18) 


-0.205(20) 


0.83(7) 


-0.37(13) 


1.16 


35.62(7) 


3.32(7) 


-0.08(8) 


-0.162(17) 


0.50(5) 


-0.07(6) 


1.23 


37.73(7) 


2.98(7) 


-0.13(6) 


-0.140(22) 


0.35(5) 


-0.09(5) 


1.36 


39.47(5) 


2.67(4) 


-0.03(3) 


-0.063(18) 


0.19(3) 


-0.01(2) 


1.50 


40.34(4) 


2.54(3) 


-0.06(2) 


-0.043(16) 


0.15(2) 


-0.03(2) 


1.65 


40.81(4) 


2.40(2) 


-0.04(2) 


-0.029(14) 


0.10(1) 


-0.02(2) 


1.81 


41.05(3) 


2.37(2) 


-0.04(2) 


-0.040(14) 


0.11(2) 


-0.01(1) 


1.98 


41.20(3) 


2.29(2) 


-0.00(1) 


-0.051(13) 


0.08(1) 


0.02(1) 



Table 3.3: Taylor expansion coefficients c7(T) and c!;f{T). 



identical strength also in x ■ This, in turn, imphes that the singular behavior will add up 
constructively in the quark number susceptibility whereas it can cancel in the isovector 
susceptibility giving rise to finite values for xi such a critical point. Even for smaller, 
non-critical values of fJ^q/T, however, the rapid rise of df'{T) for T ~ Tq is important. It 
shows that non-diagonal susceptibilities will become large at non-zero chemical potential 
in the transition region from the low to the high temperature phase, i.e. fluctuations in 
different fiavor channels, which are uncorrelated at high temperature, become strongly 
correlated in the transition region. This correlation is also reflected in the errors of the 
various expansion coefficients, which are of similar size for c^" and c!^^ but much reduced 
in the difference, xi- 

The above considerations also suggest that the electric charge susceptibihty, 

g (4Xu« ~l~ Xdd '^Xud) 

\{xi + lx?l , (3.9) 
will be singular whenever the diagonal and non-diagonal susceptibilities are singular as 
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Figure 3.4: The Taylor expansion coefficients c™ (upper row) of diagonal and c^"^ (lower 
row) of non-diagonal susceptibilities for n = 2, 4 and 6 as functions of T/Tq. 



the cancellation between the corresponding singular parts will be incomplete. We show 
the charge susceptibility in Fig. ()3.5|1 . As it is dominated by the contribution from the 
isovector susceptibility any possible singular contribution arising from Xq will be weak. It 
thus may not be too surprising that a peak does not yet show up in xc- 



3.3 Mass derivatives and chiral condensate 

The transition between low and high temperature phases of strongly interacting matter is 
expected to be closely related to chiral symmetry restoration. It is therefore also of interest 
to analyze the dependence of the chiral condensate on the quark chemical potential. We 
will do so in the framework of a Taylor expansion of the grand potential, 



rp3 



Nr\^ dlnZ 
IfJ dm IT 



ri=0 



T 



(3.10) 



with 
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Figure 3.5: The charge susceptibihty Xc/T"^ as a function of T/Tq for various Hq/T ranging 
from /iq/T = (lowest curve) rising in steps of 0.2 to /iq/T = 1, calculated from a Taylor 
series in 6*^ order. 



Here we expressed the bare lattice quark masses, ma, in units of the temperature by using 
m/T = maNT-. For n > the expansion coefficients of the chiral condensate, are directly 
related to derivatives of the expansion coefficients of the grand potential Q with respect 
to the quark mass, i.e. cf^^ = dcn/d{m/T). For n = this holds true up to a contribution 
arising from the normalization of the pressure at (T = 0, /i^ = 0). As such the coefficients 
cf^^ also provide information on the quark mass dependence of other thermodynamic 
observables like pressure, number density or susceptibilities. For instance, the change of 
the quark number susceptibility with quark mass is given by 

"^^^ - 2cf + 12cf f ^ y + . (3.12) 



dm/T ' ^ \T 



We have calculated the derivatives of c„ with respect to the quark mass for n = 0, 2 and 4. 
These derivatives are shown in Fig. ()3.6p together with the corresponding derivatives for 
the expansion coefficients of the isovector susceptibility, which have a similar temperature 
dependence, 

- 2c^''^^ + 12ci'^^ { + 0{fi^) . (3.13) 



dm/T ' " VT 

We note that the expansion coefficients c^"^ are negative for n > and T < 0.96To. The 
chiral condensate thus will drop at fixed temperature with increasing fiq/T and the chiral 
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Figure 3.6: The Taylor expansion coefficients c^"^ of the chiral condensate for n = 0, 2 
and 4 as functions of T/Tq. Also shown are the coefficients cj^'"^^ for n = 0,2 which define 
the quark mass derivatives of the isovector susceptibility in analogy to Eq. ()3.12|) . 



susceptibilities will increase in the hadronic phase with decreasing quark mass. This, 
together with the change of sign in cf^ at T ~ Tq, will shift the transition point at non- 
zero fiq/T to lower temperatures. In Fig. ()3.7|) we show the chiral condensate and the 
related chiral susceptibilityj obtained from a Taylor expansion up to and including C(/iq), 



2^2 



trM" 



(3.14) 



Obviously x^iv develops a much more pronounced peak for Hq/T > than at vanishing 
chemical potential which, moreover, is shifted to smaller temperatures. However, as 
will become clear from the discussion in the next section the peaks found in and 
also in other susceptibilities have to be analyzed and interpreted carefully. They reflect 
the abrupt transition from the hadronic regime to the high temperature phase in which 
fluctuations of the chiral condensate are suppressed, but do not signal the presence of a 2"'^ 
order phase transition unambiguously. The rapid rise of susceptibilities in the hadronic 
phase is strongly correlated to the increase in the pressure and is also present in a hadron 
gas which does not show any singular behavior at the transition temperature. 

Also the expansion of the chiral condensate and related observables are compatible 
with the HRG model. In fact, a comparison of the temperature dependence of c^"^ shown 
in Fig. ()3.(')j) with that of the expansion coefficients c„ of the grand potential shown in 

J The chiral susceptibihty introduced here is not the complete derivative of the chiral condensate with 
respect to the quark mass. As frequently done also at — we define the chiral susceptibility by ignoring 
a contribution from the connected part which would arise in the derivative d(ip'>p) /dm. Nonetheless 
seems to capture the leading singular behavior that should show up at a 2'"^ order critical point [30] . 
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Figure 3.7: The chiral condensate {'ipip) (left) and chiral susceptibility Xij-ip (right) as a 
function of T/Tq for /ig/T = 0, 0.4 and 0.8. The chiral condensate drops with increasing 
Hq/T and the peak in becomes more pronounced. 



Fig. ()3.1|) suggests a strong similarity between —cf^^ and Cn+2- On the other hand, for 
T < To the ratio cf^ /c2^ agrees within errors with the ratio C4/C2 shown in Fig. ()4.H1 . All 
this is consistent with the HRG model where the quark mass (or spectrum) dependence 
only enters through the functions F{T) and G{T) and does not modify the dependence 
on fig/T. 



4 Radius of convergence and the hadron resonance 
gas 



So far we have not discussed the range of validity of the Taylor expansion. In general 
the Taylor series will only converge for fiq/T < p (or fiq/T < p) where the radius of 
convergence, p, is determined by the zero of Z(T, pg, pg) closest to the origin of the 
complex fig plane. If this zero happens to lie on the real axis the radius of convergence 
coincides with a critical point of the QCD partition function. A sufficient condition for 
this is that all expansion coefficients are positive [31]. Apparently this is the case for all 
coefficients c„(T) with T/Tq < 0.96 that have been calculated so far by us*^. Above Tq, 
however, we find from the calculation of Cq{T) that the expansion coefficients do not stay 

'^In [20] it is reported that Cs{T) is negative for T < 0.95To, however, the statistical significance of 
this result unfortunately is not given. 
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Figure 4.1: (a) The ratios C4/C2 and C6/C4, and (b) C4/C2 and cl/c{ as well as (c) /cq^ 
and cl^ /c2^ as functions of T. Horizontal lines indicate the HRG prediction for T < Tq, 
and the SB prediction for T > Tq. Note the difference in vertical scale between the plots. 



strictly positive. This is in accordance with our expectation to find a chiral critical point 
at some temperature T <Tq. 

The radius of convergence of the Taylor series for Q{T, fig, fig) can be estimated by 
inspecting ratios of subsequent expansion coefficients, 



p = lim p2n = lim 



C2r, 



C2n+2 



(4.1) 



where the square root arises because the Taylor expansion of the grand potential VL is an 
even series in Pq/T. The ratios C2n+2/c2n are shown in Fig. ()4.H) together with ratios of 
the expansion coefficients of the isovector susceptibility and cf^^ of the chiral conden- 
sate. It is obvious that these ratios rapidly change across Tq and approach the value of 
corresponding ratios obtained in the high temperature ideal gas limit. Another remark- 
able feature, however, is that below Tq the ratios involving expansion coefficients of the 
/ig-dependent parts of fi, xi and Xvii/> almost temperature independent. In fact, these 
ratios are consistent with the corresponding ratios deduced from the grand potential of a 
hadron resonance gas (Eq. (j2.23|) ). i.e. C4/C2 = c^'^/o^'^ = 3/4 and 05/04 = c^/cl = 3/10. 
In ratios that contain the lowest order expansion coefficients, i.e. cq, Cq^ and Cg, the 
spectrum dependence does not cancel because the lowest order expansion coefficients also 
depend on the meson sector which is not the case for higher order coefficients. These ra- 
tios thus show a significant temperature dependence as can be seen for c{/c2 and cp' /cq^ 
shown in Fig. ()4.1|) . 

Similar information is contained in the ratios of physical observables, e.g. the quark 
number density or pressure over the quark number susceptibility, introduced in Eq. ()2.20p . 
These ratios are shown in Fig. ()4.2|1 . Here UgjT^ and XqjT'^ have been calculated using 
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Figure 4.2: The baryonic part of the pressure divided by the quark number susceptibihty, 
Ap/xq (left), and the normahzed derivative of pressure with respect to quark number 
density, Ug/xq (right), as a function of fiq/T for various T/Tq. Horizontal lines show the 
infinite temperature ideal gas values and the HRG model prediction (T < Tq) (solid lines) 
and expanded to 6*/i-order in Hq/T (dashed lines). The difference is visible only for Ap/xq 
at fijT = 0.8. 



Eqs. fl3.3|) and ()3.4|) up to 0(/ig). The ratio, riq/xq = (dp/dfiq) /{duq/dfiq) = dp/driq is 



related to the isothermal compressibility, kt = Xql f^q which diverges at a 2 order phase 
transition point, i.e. at a point at which dp/dnq = 0, the number of particles is unstable 
under small changes in the pressure (mechanical instability) and large density fluctuations 
occur. This instability leads to a divergence in the quark number susceptibility [21]. A 
2nd QTi-(^Qj- phase transition is thus expected to be signalled by a zero in both ratios shown 
in Fig. ()4.2|1 . On the other hand, for /i/ = these ratios are expected to be constant in 
an ideal quark-gluon plasma as well as in a hadron resonance gas. 



^SB 



tanh [-^] . (4.2) 



The corresponding values are indicated in Fig. ()4.2|) by horizontal lines. 

As far as the determination of a possible 2"'^ order critical point at non-vanishing 
quark chemical potential (chiral critical point) is concerned Fig. (j4.ip and Fig. (j4.2j) 
contain identical information. For T < 0.96To bulk thermodynamic observables agree 
with predictions based on an HRG model, which in itself does not show any critical 
behavior as function of fiq/T at fixed T. In particular, there is no hint for a dip in 
^p/Xq or ^g/Xq which could signal the presence of a second order transition point. The 
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same observation, albeit with larger statistical errors, holds for ratios involving the chiral 
susceptibility x^^v Nonetheless, all these quantities change rapidly in the transition from 
the low temperature to the high temperature regime and, moreover, at T = Tq the 6*^ 
order expansion coefficients clearly cannot be described within the HRG model. Due to the 
good agreement with the HRG model and its Taylor expansion at lower temperature we 
cannot, however, present an upper limit for the radius of convergence below Tq; the ratios 
shown in Fig. (j4.1|) suggest that a lower limit is given by (/Xq/T)c>l. Also from the analysis 
of the temperature dependence of bulk thermodynamic observables we get, at present, no 
unambiguous evidence for the existence of a phase transition. At present, therefore, we 
cannot rule out that in the temperature range covered by our analysis (T>O.STq) the 
transition to the high temperature phase is a rapid crossover transition rather than a 
phase transition. This situation then would be similar to that at /ig = 0. In order to 
exclude this possibility we would need, in the future: to consider even higher orders in 
the Taylor series; to scan in more detail the small temperature interval [0.95To, Tq]; and to 
explore systematically the quark mass and volume dependence of our results. These issues 
are partially addressed already in the next section where we discuss the use of reweighting 
techniques to calculate some thermodynamic observables and compare results obtained 
within this approach with results from the Taylor expansion. 

The good agreement found here for different ratios of Taylor expansion coefficients 
calculated on the lattice and within the HRG model suggests that we may use this infor- 
mation for a more detailed analysis of the composition of hadronic matter at temperatures 
below Tq. In Ref. [15] also the temperature dependence of thermodynamic observables like 
the pressure or the quark number susceptibility have been compared to the HRG model.' 
In order to do so the hadron spectrum has been adjusted to the conditions realized in 
the lattice calculations, i.e. all masses have been shifted to larger values as the lattice 
calculations have been performed with unphysically large quark masses. This approach 
can also be turned around. The HRG model for T < Tq can be used as an ansatz to 
determine the contributions of the mesonic and baryonic parts of the spectrum without 
making assumptions on the distortion of the spectrum due to the unphysical quark mass 
values. 

As outlined in section 2, within the Boltzmann approximation the HRG model yields 
a simple dependence of the pressure on the quark chemical potential. The relation given 
in Eq. ()2.15p can easily be extended to also include a non-vanishing isovector chemical 
potential. Neglecting the mass difference among isospin partners the pressure can be 
written as. 



PJZl^ ^ G« (T) + G(3)(T) 1(2 cosh 



For these observables a good functional agreement between lattice data and the leading dependent 
term in the HRG model has also been noted in [14] within the imaginary chemical potential approach. 
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Figure 4.3: Contribution of different hadronic channels to the total pressure p/T^ obtained 
by using the HRG ansatz. Shown are results for (a) fiq/T = 0, (b) fiq/T = 0.4 and (c) 

fXg/T = 0.8. 



+F(^)(T) cosh (^) cosh (^) 
+F(^)(T)icosh(^) [cosh(^ 



cosh (j;-^ + cosh (^^^ 



(4.3) 



where G^-^\ G'^^\ F^"^^ and F*^^) are the contributions to the pressure at = /x/ = arising 
from isosinglet mesons (//,..., [Bi = 0, = 0]), isotriplet mesons (vr, . . . , [Bi = 0, I^i = 
{0, ±1}]), isodoublet baryons {n,p, . . . , [B^ = ±1, I^i = {±1/2}]) and isoquartet baryons 
(A,..., [Bi = ±l,/3i = {±1/2, ±3/2}]), respectively. These functions contain all the 
information on the hadron spectrum in different quantum number channels. Performing 
the Taylor expansion of the pressure as well as quark number and isovector suscepti- 
bilities allows to relate these functions to combinations of the various Taylor expansion 
coefficients. This way one finds 



G'(3)(T) 
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F(2)(T) 
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(4.4) 



The various contributions to the pressure are shown in Fig. (j4.3p a for /Xg = 0. With 
increasing quark chemical potential the relative weight of hadrons in different quantum 
number channels changes. As expected the baryonic component becomes more important 
with increasing fig (Fig. ()4.3p b and c). We find that for fj,q/T>0.6 the baryonic sector 
gives the dominant contribution. 
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Figure 5.1: Contour plot of the variance of the phase of the quark determinant, a{6) 
calculated for 6**^^^ in the {T/Tq, fXq/T) plane. Contour lines for a{9^^^) are given in steps 
of 7r/4 ranging from ir/A (lowest curve) to 27r. 



5 Re weighting Approach 



An alternative to a strict Taylor expansion of thermodynamic observables in terms of 
^q/T is the reweighting approach. Here the dependence of the grand potential on the 
quark chemical potential is included in the calculation of observables, X, by shifting the 
/Xq-dependent piece of the QCD action into the calculation of expectation values rather 
than taking it into account in the statistical weights used for the generation of gauge field 
configurations. This reweighting approach has been used to analyze the thermodynamics 
of QCD at non-zero chemical potential [10,12]. Within this approach thermodynamic 
observables X{P,fj,) are estimated via the expression 



/Y\ _ {X ei e ViPofl) , . 



where AlndetM = lndetM(/i) - IndetM(O) and ASg = Sg{p) - Sg{po) is the difference 
of the gluonic part of the QCD action. The expectation values on the RHS of Eq. 1)5.11) are 
obtained in simulations at {j3o,0). In [3] we implemented a version of Eq. ()5.H) in which 
the reweighting factor A In detM as well as the operator X itself have been replaced by 
a Taylor series about n = 0. The advantage over an exact evaluation of detM [2] clearly 
is that the required expressions are calculable with relatively little computational effort 
even on large lattices. In our initial study we performed the expansion consistently up 
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Figure 5.2: Susceptibilities Xg/T"^ (left) and Xi/T'^ (right) for various fiq/T ranging from 
/ig/T = (lowest curve) rising in steps of 0.2 to fiq/T = 1. Results are obtained by 
a combination of reweigliting (solid lines and data points) and from a 6*''-order Taylor 
expansion (dashed lines). 



to and including (9(yU^). Here we extend this analysis by expanding IndetM up to and 
including terms of (9(/i®). Unlike the direct evaluation of thermodynamic observables in 
terms of a Taylor expansion up to a certain order the reweighting approach with a Taylor 
expanded weight factor also includes effects of higher orders in fiq/T which are partially 
resummed in the exponentiated observables exp{-^A In detM}. 



The effectiveness of any reweighting approach strongly depends on the overlap between 
the ensemble simulated at (/^c/^o = 0) and that corresponding to the true equilibrium 
state at (/?, fi) which one wants to analyze. This can be judged by inspecting the average 
phase factor of the complex valued quark determinant, (e*^)(/3Q^o)) where the quark deter- 
minant is written as detM = |detM|e*^. Reweighting loses its reliability once (e*^) ^ 1 as 
both expectation values appearing in the numerator and denominator of Eq. ()5.H1 then 
become difficult to control [32]. In our approach we estimate the phase factor via the 



variance of the phase 9, a(6) = y {O"^) — {Oy, where we approximate the phase by its 
Taylor expansion up to (9(yU^"~^), 

g(n)_^fl^f c^^^-MndetM 

As discussed and shown in Fig. 6 of [32], the value of /ig/T for which the standard deviation 
of 9^""^ exceeds 7r/2 is a reasonable criterion for judging the applicability of reweighting 
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in our simulated systems. In Fig. ()5.1|) we show contour lines for the variance of 9^^\ 
All contour lines drop dramatically in the vicinity of Tq; the contour corresponding to 
a{9^^^) = 7r/2 yields /ig/T ^ 1.5 at T ~ 1.2To and reaches a minimum value of about 0.6 
at T ~ 0.9To. Reweighting is thus much easier to control in the high temperature phase 
than in the hadronic phase. 

In Fig. ()5.2|) we plot quark number and isovector susceptibilities for various fiq/T cal- 
culated by using reweighting where possible. The operators needed for the susceptibilities 
are calculated with an error of 0{fi^), the reweighting factor A IndetM with error 0{fi^). 
Dashed lines show a comparison with the results obtained from direct Taylor expansion 
of the susceptibilities up to and including (see Fig. ()3.3p ). 

As anticipated above, for T<Tq reweighting becomes difficult to control for yUg/T > 0.6. 
However, where reweighting appears to be statistically under control, it agrees well with 
the direct Taylor expansion and even shows similar features in the statistical errors; i.e. 
the signal is much noisier in the vicinity of Tq for Xq than for xi- 



6 Conclusions 

We have extended our analysis of the thermodynamics of two flavor QCD at non-zero 
quark chemical potential to the 6*^ order in a Taylor expansion around fiq/T = 0. We 
find clear evidence for a rapid transition from a low temperature hadronic phase to a high 
temperature quark-gluon plasma phase which is signalled by large fiuctuations in the 
quark number density and the chiral condensate. The transition temperature shifts to 
smaller values with increasing quark chemical potential. Above Tq the Taylor expansion 
coefficients and bulk thermodynamic observables agree with qualitative features of the 
perturbative high temperature expansion and approach the ideal gas limit to within ~ 20% 
for T>1.5Tq. Thermodynamics in the low temperature phase agrees well with predictions 
based on a hadron resonance gas for temperatures T<0.96To and fiq/T<l. 

From the analysis of bulk thermodynamic observables alone we cannot provide strong 
evidence for the existence of a 2"^ order phase transition point in the QCD phase diagram. 
At present we cannot rule out the transition being a rapid crossover in the entire parameter 
space covered by our analysis. In particular, we have shown that large fluctuations in the 
quark number and the chiral condensate are consistent with expectations based on a 
hadron resonance gas. The current estimates on the radius of convergence of the Taylor 
expansion favor a critical value of the quark chemical potential close to fiq ~ Tq. The good 
agreement of the expansion coefficients with those of a hadron gas, however, prohibit any 
firm conclusion on the location and even on the existence of the chiral critical point. 

Likewise, we cannot rule out that a 2"^^ order transition occurs at temperatures closer 
to To than the largest value in the hadronic phase, T = 0.96To which we have analyzed 
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here. In order to improve on the current analysis it would be important to perform 
calculations at smaller quark masses in a narrower temperature interval around Tq. An 
improvement over the current statistical errors on the G*'^ order coefficient as well as the 
analysis of higher order expansion coefficients with high statistics is needed. 
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A Appendix: Taylor expansion coefficients 

In this appendix, we derive some equations which are used in the calculation of the various 
thermodynamic quantities and expansion coefficients of the Taylor series presented in this 
study. The partition function is given by 

Z = yPC/(detM)"f/V*« , (A.l) 
with U e SU{3). The expectation value of a physical quantity, (C) is then obtained as 

(C») = ^yP[/C>(detM)"f/V^« , (A.2) 
and its derivatives with respect to quark chemical potential and quark mass are given by, 

d{0) ^ ldO\ nf / / 9aiidotM) \ _ / 9(hidetM) \\ 
dm \dml A \\ dm / \ dm /) 

Here we use m as the dimensionless quark mass value instead of ma, and also ^ = yU^a 
for the dimensionless quark chemical potential. The temperature is T = (AT^a)"^ and the 
volume is V = (N^a)^. Moreover, we introduce for simplification, 

_ rif d'^trM-^ _ m 9"+^ In det M ^ _ rif 9" In det M 
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All Taylor expansion coefficients used in this paper can be expressed in terms of expec- 
tation values of certain combinations of C„ and X>„. The required derivatives of IndetM 
and trM~^ are explicitly given in the following. 



Derivatives of In det M: 



91ndetM / , ,dM\ 

92 IndetM / ^ ,d'^M\ / ^ ,dM ,dM\ 

IndetM _ _ ( M-^^M-^^-^\ 

+2tr M-^^M-i^M-i^ , A.8 
y a/x Of! on J 

-6tr M-'--M-'^M-'^M-'—- , A.9 
y a/x o/x a/i c/yu J 

+30tr M-'--M-'—^M-'—^ 

.r. ( i-^M^ i^M^ i^M^ ,a2M\ 
-60tr M-i^M-i^M-i^M-i— ^ 

f.. i-^M,^ i«9M,^ i«9M,^ i^M,^ T(9M\ 
+24tr M-^^M-^^M-^^M-^^M-^—- , (A.IO 
y a/x a/i oil oil oil J 

a^lndetM / , i(96M\ ^ / , i'9M,^ ,d^M\ 

= " r i^j - i)^" 

y dfj, oil oii^ J \ oil 011^ oii'^ J 

+60tr M-^— -M-I--M-I— - + 30tr M"!— -M"^— -M"^— - 
\ a//^ a// a//'^ J \ on^ on^ on^ J 
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( 1^^,, 1^^,, 1^^,. ,d^M\ 
-120tr M-^--M-^--M-^--M-^-— 

\ on on on on j 

( ^^M i«9M_ i^^M^ ,d'^M\ 
-180tr ^M-i— ^ 

y on on on on J 

-90tr M-'--M-'-—M-'--M-'-— 

\ On on on on J 

+360tr M-^^M-^^M^^^M-^^M-^—- 

\ On on (yf^ ^^i^ (^f^ J 

i9M^^ i<9M,^ i<9M_ i9M_ ,dM\ , 

-120tr M'^^M-^^M-^^M-^^M-^^M-^—- . A.ll) 

y a/i on on on on on J 



Derivatives of trM ^: 



\ On on I \ on on on ) 

y a//^ a//^ y \ on on ) 
-12tr M-^—-M-^^M-^^M-^ 

\ on On On J 

-12tr M-'--M-'—^M-'--l 

\ on on on 

( i9M^^ i<9M,^ .d^M^ 
-12tr M-^^M-^^M-'—^I 

\ on On on j 

+24tr M-^^—M-^^—M'^^—M-^^—M-^ 

\ On on on on 



Having defined the explicit representation of C„ and we now can proceed to define the 
expansion coefficients for various thermodynamic quantities discussed in this paper. 
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Pressure (p) : The pressure is obtained from the logarithm of the QCD partition 
function. Its expansion is defined in Eqs. ()3.1|) and ()3.2|) . The leading expansion coefficient 
Co is given by the pressure calculated at /i^ = 0. All higher order expansion coefficients 
are given in terms of derivatives of In Z. 



^4 yrj^s nyrj. 



with 

1 9" In Z 



(A.17) 

At=0 



To generate the expansion we first consider derivatives of hiZ for /i 7^ 0. For the first 
derivative we find 

-^^A, = {V,) . (A. 18) 

Higher order derivatives are generated using the relation 

^ = A+i-AA , (A. 19) 
oil 

where An is defined as 

A - (exp{-Po}^:^^|f^) . (A.20) 

With this we can generate higher order derivatives of In Z iteratively using 

a"+MnZ/9/i"+^ = . (A.21) 

Explicitly we find from Eq. ()A.20jl 

A2 = (P2) + (Vf) , (A.22) 



^3 = (^^3) + 3 {V^V,) + (Vl) , (A.23) 

A = (V,) + 4 {V;V^) + 3 (vl) + 6 {V2VI) + (V{) , (A.24) 
A = {^5) + 5 (V^V,) + 10 {VsV2) + 10 (VsVl) + 15 (vjVi) + 10 (V2VI) 

+ {VI), (A.25) 
A = (^^e) + 6 + 15 + 10 (vl) + 15 (VJ^I) + 60 {V^V2V^) 

+15 (l)^) + 20 (^3^^?) + 45 (1)221)?) + 15 (r-s^^f ) + . (A.26) 
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From Eq. ()A.21|) we then obtain through repeated apphcation of Eq. ()A.19|) . 
d\nZ 



(9/i 

aMnZ 
d^nZ 



Ai, 
A2 - Al 



(A.27) 
(A.28) 

= ^3 - 3^2 A + ^Al (A. 29) 

= A4- AA3A1 - 3Al + 12^2^? - QAt, (A.30) 
= ^5 - 5AiAi - IOA3A2 + 20^3^? + 30^2^1 - 60^2^? + 2AAl (A.31) 
= As- 6A5A1 - 15A^2 - 10^3 + SOAiAj 

+ 120^3^^! + 30^^ - 120^3^? - 270^^^? + 360^^1 - 120^? . (A.32) 



These relations simphfy considerably for /i = as all odd expectation values vanish, i.e. 
An = for n odd. In fact, (9"'(ln det M) / dfi^ is strictly real for n even and pure imaginary 
for n odd. Using this property, the odd derivatives of the pressure vanish and also the even 
derivatives become rather simple. This defines the expansion coefficients c„ introduced in 
Eqs. dSH) and (Q, 



C2 



C4 = 



C6 



1 d\p/T^) 


tlq=Q 


_IN 


1 d\p/T^) 




1 1 


- 4! ^{^^JTY 




4! N^N^ 


1 d^{p/T^) 




1 1 


6! di^ijTf 


fiq=0 


6! N^N^ 



-{Ai-3Al) , 
■{Ae - I5A4A2 + 3OAI) 



(A.33) 



Here all expectation values An are now meant to be evaluated at /i = 0. 



Isovector susceptibility (xi)'- While the Taylor expansion of the quark number sus- 
ceptibility is easily obtained from that of the pressure we need to introduce the expansion 
of the isovector susceptibility. This has been done in Eq. ()3.5|) . More explicitly the 
isovector susceptibility is given by 




dHnZ dHnZ 



la^QndetMj 
4 dj?f 




S^lnZ d'^lxiZ' 



1 a(ln det M„) 
4 dfiu 



1 (9(ln det Md 
4' 
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'l g(lndetM^) _ i a(lndetMrf) \^ 
4 djiu 4 djid I 



(A.34) 



Throughout this paper we have considered the case of flu = fid = fJ'qd = fJ^- The 
isovector chemical potential /// has been set equal to zero after appropriate derivatives 
have been taken. In the Taylor expansions of Xh which is performed at /x = in 
terms of fJ^q/T, the derivatives with respect to /i„ and fid then become identical, i.e. 
[d-{\ndetMu^d))/dK^d)]iMd)) ^ ^ (In det M)/9/i"](/i). 



The calculation of the isovector susceptibility thus reduces to the calculation of 
Xi Nr /52(lndetM)\ 



T2 2N^ \ dfi^ 



(A.35) 



To define the expansion of the isovector susceptibility, xi at fixed /i/ = around fiq = 
we set up an iterative scheme similar to that introduced for the pressure. We introduce 
the additional kernel V2 in Eq. ()A.20|) to generate Bn for arbitrary /ig 7^ 0, 

Bn+2 = ( exp{-Po} ) ■ (A.36) 

This yields. 



B2 = {V2) , (A.37) 

Bs = m + {V^V^) , (A.38) 

= {V^) + 2 {V^V^) + {Vf) + {V2VI) , (A.39) 

B, = {V,) + 3 (V^V,) + 4 {V,V2) + 3 (V.VI) + 3 (vlV,) + (v^vf) (A.40) 

= m + 4 {V^V^) + 7 (P4^^2) + 4 (1)2) + 6 {vj)l) + 16 (I?3^)2^)i) 

+3 (Vl) + 4 (1)3^^?) + 6 Ivlvf) + (P2^)^) . (A.41) 



The derivative of Bn with respect to the chemical potential satisfies, 
dBn 



Bn+i-BnAi . (A.42) 



Starting with {N^/Nl)B2 = Xi/T^, we obtain the equations for d''{xi/T^)/d{i2g/T)'' 
iteratively, and then use again the CP symmetry for fi = 0, i.e. An and Bn are zero for n 
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odd. Therefore the odd derivatives in the expansion of XijT"^ at /i = are zero and the 
even derivatives define the expansion coefficients used in Eq. fj3.5|) . 



1^ 

2T2 



^lq=0 



IN 



1 d'iXi/T' 



4! difijrr 

1 d\xi/T') 



Hq=0 



6! difijTY 



1 1 
1 1 



^,=0 6! iV3iV3 



{B,-B2A2), 
{Be - QB4A2 



B2AA + QB2AI) 



(A.43) 



Chiral condensate {(xfyt/j)) and disconnected chiral susceptibility (x^^^): We 

also discuss the Taylor expansion of the chiral condensate and the related chiral suscep- 
tibility, 



X'tpip 
2^2 



_^ / glndetM \ _ 

dm ~ iV|T ^ ~ ^ 



N^dlnZ _ N^i 



n clet M \ r^^rii / ^ 
dm / 4 \ 



A 



/(trM-i)^\ - (trM-i) 



(A.44) 
(A.45) 



The iterative scheme is similar to that introduced for the isovector susceptibility but with 
the generating kernel P2 in Eq- ()A.36|) replaced by Cq. This yields 



•^0 — (^^o) ; 

^1 = (Ci) + (CoPi) , 

^2 = (C2) + 2 (CiDi) + (C0P2) + {CoVl] 

^3 = (C3) + 3 (C2P1) + 3 (C1P2) + {CoV^) + 3 ICiVl) + 3 {C0V2V1) 

^4 



(A.46) 
(A.47) 
(A.48) 

(A.49) 

(C4) + 4 (CsDi) + 6 {C2V2) + 4 (C1P3) + (CoI?4) + 6 (C2P1') + 12 (CiP2^^i) 
+4 (CoPs^^i) + 3 IC^vf) + 4 (CiP?) + 6 lCoV2Vf) + (CqP^) . (A.50) 



These expectation values again satisfy 
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(A.51) 



This gives the expansion coefficients c^"^ for the chiral condensate defined in Eqs. ()3.10p 
and (jTTT|l . 



Co 



C2 



C4 



T3 



fj.q=0 



1 92 /T3 



Atq=0 



2 9(/x,/T)2 
1 (^^) /T3 



1 1 d'^U^)/T^ 



2N3 



1 1 d^Ui^/T"^ 



flq=0 



Aim 



4! 9(^,/T)4 



M=0 



2 ' 



(A.52) 



(A.53) 



(A.54) 



where we used again that An and JF„ are zero for n odd ant /i = 0. Hence the odd 
derivatives in the expansion vanish. Note that these expansion coefficients also control 
the quark mass dependence of the quark number susceptibility given in Eq. (j3.12j) . The 
first derivative of the isovector susceptibility with respect to the quark mass, 



dm/T 



1 

iv! 



'uf 93 (In det M) 
^ 4 dfi^dm 



4/ (9^2 ^ ^ 



/ n, {In det M) \/ni, 
\ 4 9/i2 / \ 4 



(A.55) 



for which we have introduced the Taylor expansion in Eq. ()3.13|) requires the calculation 
of further expectation values In, 



T3 = (C3) + (C2P1) + (CiV^) + (C0P3) + {C0V2V,) , 

J4 = (C4) + 2 (CsDi) + 2 (C1P3) + 2 (C2P2) + (C2l^?) + 2 (CiP2^^i) 



+ (C0P4) + 2 (CoI?3l?i) + (Co2^: 



C0V2VI 



(A.56) 
(A.57) 

(A.58) 



With these and the coefficients Bn and JF„ one finds for the expansion coefficients, 



'-2 



C4 



1 9(x//T2) 



2 dm/T 



/.=o 2 A^ 



fj.=0 



4! d{^ijTYdm/T 

\\WN^^^^ - X2A - ^4-^0 - i32 J-2 + 2^2-^0 A) 



(A.59) 
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Finally we present the expansion of x^^) which is generated in analogy to the isovector 
susceptibility by using Cq as a kernel in Eq. (|X20|l . We find 



6^0 = _ 

e?i=2(CiCo) + (Co'Pi), 

= 2 (C2C0) + 2 {Cl) + 4 (CiCoPi) + {CIV^) + (Co^P?^ , 

= 2 (C3C0) + 6 (C2C1) + 6 (CsCoI^i) + 6 (C^Di) + 6 {C^C^V^) + (Cq'I^s) 

+6 (CiCoP?) + 3 {CIV^V^) + (C^P?^ , 
= 2 (C4C0) + 8 (C3C1) + 6 (Ci) + 8 (CsCoVi) + 24 (Csdl^i) 

+12 (CaCoDs) + 12 (C^Ps) + 8 (CiCqI^s) + (C^V,) + 12 (CaCoD?) 

+12 (C^Vl) + 24 {CiCoV^V^) + 4 (Cq^Ps^^i) + 3 (C^qVI) + 8 (CiCqI)?' 

+6(Co^I),I)?) + (Co^Pt). 



(A.60) 
(A.61) 
(A.62) 

(A.63) 



(A.64) 



These expectation values satisfy, 



dGn 
9/i 



(A.65) 



The expansion of the disconnected chiral susceptibility defined in Eq. (j3.14|) is then given 

by, 



^2 



2^2 



(6^0 



■2^ 

; ' 



2 a(/i,/T)2 
1 d%^/T' 



[Q2 - Q0A2 - 2 - Jo^2) ^0 



4! 9(/i,/T)4 ^=0 
, \ 3 [e?4 - 6e?2^2 - Q^Ai + 66^0^2 - 6 (^2 - -^0^2)' 

a T 

-2 (^4 - - :FqAa + 6 ^0^2) -^0 



(A.66) 
(A.67) 



(A.68) 
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